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Abstract 

The mechanism of DC-Electric-Field-Induced Second-Harmonic (EFISH) 
generation at weakly nonlinear buried Si(001)-SiO2 interfaces is studied ex- 
perimentally in planar Si(001)-SiO2-Cr MOS structures by optical second- 
harmonic generation (SHG) spectroscopy with a tunable Ti:sapphire femtosec- 
ond laser. The spectral dependence of the EFISH contribution near the direct 
two-photon E\ transition of silicon is extracted. A systematic phenomeno- 
logical model of the EFISH phenomenon, including a detailed description of 
the space charge region (SCR) at the semiconductor-dielectric interface in 
accumulation, depletion, and inversion regimes, has been developed. The in- 
fluence of surface quantization effects, interface states, charge traps in the 
oxide layer, doping concentration and oxide thickness on nonlocal screening 
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of the DC-electric field and on breaking of inversion symmetry in the SCR is 
considered. The model describes EFISH generation in the SCR using a Green 
function formalism which takes into account all retardation and absorption 
effects of the fundamental and second harmonic (SH) waves, optical interfer- 
ence between field-dependent and field-independent contributions to the SH 
field and multiple reflection interference in the Si02 layer. Good agreement 
between the phenomenological model and our recent and new EFISH spec- 
troscopic results is demonstrated. Finally, low-frequency electromodulated 
EFISH is demonstrated as a useful differential spectroscopic technique for 
studies of the Si-Si02 interface in silicon-based MOS structures. 

PACS numbers: 42.65.Ky, 73.65. Qv, 68.35.-p 
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I. INTRODUCTION 



Optical Second Harmonic Generation (SHG) has been one of the most intensively studied 
phenomena in surface and interface optics [|l]-|3|] for the last decade. The interest in SHG 
stems from its unique sensitivity to the structural and electronic properties of surfaces and 
interfaces of centrosymmetric media. This unusually high surface/interface-sensitivity comes 
about because, in the electric dipole approximation, SHG is forbidden in the bulk of materials 
with inversion symmetry ffij , but allowed at interfaces, where inversion symmetry is broken 
by the discontinuity of crystalline structure. Related nonlinear sources of SHG are localized 
in a thin (several nanometers thick) surface or interface layer. In semiconductors, inversion 
symmetry is also broken by the DC-electric Field (DCF) in the subsurface Space Charge 
Region (SCR), which is created by initial band bending and/or external bias application. 
The lack of inversion symmetry in the SCR results in DC- Electric-Field Induced Second- 
Harmonic (EFISH) generation, which manifests itself through electromodulation of the SHG 
intensity. Thus, all important properties of surfaces, buried interfaces and subsurface layers 
- their charge electronic surface state density [PHTT|, roughness (morphology) fl2| , |T3|l , 

adsorption (adatom and admolecule surface density) |14|-- |T7j1 , initial band bending |18H20|, 
etc. - can, in principle, be determined by means of the SHG probe. 

The technological importance of Si(001)-SiO2 interfaces stems from their ubiquitous 
presence in Metal-Oxide-Semiconductor (MOS) structures and MOS Field Effect Transis- 
tors (MOSFET). EFISH generation provides a promising noninvasive, in situ technique 
for characterizing interfacial imperfections and charge defects at the Si(001)-SiO2 interface. 
Moreover, the relative simplicity of the description of the SHG response from the Si(OOl) 
face, originating from the small number of tensor components of the interface quadratic 
susceptibility and the rotationally isotropic interfacial SHG response, makes Si(001)-SiO2 
interface among the most important for investigation of fundamental aspects of the EFISH 
phenomenon. 



The 1967 discovery of EFISH generation by Bloembergen and co-workers ETI] at Si- 
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electrolyte and Ag-electrolyte interfaces in electrochemical cells remained largely unnoticed 
for a number of years. The 1981 discovery of surface-enhanced SHG by Shen and co- 
workers p2[ rejuvenated interest in this effect. Surface-enhanced EFISH generation at a 



silver-electrolyte interface was observed shortly afterward ||23|| . Since 1984 EFISH has been 
systematically studied at Si(lll)- electrolyte interfaces p4]-f2^], and to a lesser extent at 
other semiconductor-electrolyte interfaces: Cd 3 P 2 (HI) H, CdIn 2 S 4 (111) |3(J, GaN (001) 



3i| , TiO"2 [p^l . These studies revealed that the strength of the DC-electric field which 
could be applied electrochemically was limited by interface electrochemical reactions, such 
as oxidation of a silicon surface at anodic potential. To circumvent this restriction, EFISH 
generation studies were extended to Si-Si0 2 MOS structures with bias applied by a ring 



metal |33| or semitransparent Cr fll~9| , |34|j gate electrode, and to GaAs-based MOS structures 



A simple phenomenological model of EFISH based on the "interface field approximation" 
- which assumes linear dependence of the DC-field-induced nonlinear polarization on inter- 
face DCF strength and yields quadratic dependence of EFISH intensity on bias voltage - was 
developed for the Si-Si02-electrolyte interface in Refs. [25,26]. Since clear deviations from 
a quadratic bias dependence were observed p3,B3], this model was improved by taking into 



account the nonlinear interference of DC-field induced and field-independent contributions 



to the nonlinear quadratic polarization as well as retardation and absorption effects pi 
Further improvement resulted from considering the spatial inhomogeneity of the DCF and 
the DC-electric-field- induced contribution to the nonlinear polarization ||34|| . These effects 
were later analyzed with a Green-function formalism p6| , |37| . At present, the most com- 
prehensive description of the EFISH phenomenon is presented in Ref. [37]. However, this 
analysis remains incomplete on three points. First, it is restricted to the depletion regime of 
the SCR, whereas experimentally applied biases have included accumulation and inversion 
regimes. Moreover, as we demonstrate in this paper, the transition from depletion to inver- 
sion to accumulation drastically changes the EFISH response. Second, surface quantization 
effects originating from strong field localization in inversion and accumulation regimes, as 



well as the role of interface states, should be taken into account. Third, multiple reflection 
interference in the Si02 layer, which significantly affects the SHG intensity from Si-SiC>2 



structures |p8|-[40|1, was neglected. 

In this paper we present a comprehensive phenomenological model of EFISH generation 
supported by experimental spectroscopic studies of p- and n-type Si(001)-SiO 2 -Cr MOS 
structures. The key features of our model are: 1) a detailed electrophysical model of the 
SCR in the accumulation and inversion regimes, which takes into account interface states 
and oxide charge traps and their effect on the spatial DCF distribution in the SCR; 2) a 
rigorous nonlinear optical model of EFISH in the SCR, based on a Green-function formalism, 
which takes into account all retardation effects, absorption of the fundamental and SH 
radiation, multiple reflection interference of both the fundamental and SH waves in the 
oxide and optical interference between field- dependent and field-independent contributions 
to the quadratic nonlinear polarization. The key feature of our experiments is comprehensive 
observation of the dependence of SHG on numerous parameters, including applied bias, 
azimuthal sample rotation, wavelength near the direct two-photon E\ transition, doping 
concentration, and oxide thickness. These combined dependences allow us to deconvolve 
the EFISH contribution fully from field-independent contributions. The non-quadratic bias 
dependence of the EFISH intensity predicted in Refs. [32,37], and its variation with doping 
concentration, oxide thickness, interfacial state density, and wavelength, is observed and 
analyzed in detail. 

II. THEORETICAL BACKGROUND 

A. Quadratic optical response of the Si-Si02 system 

In the presence of a DCF the nonlinear polarization of a centrosymmetric semiconductor 
at the second-harmonic (SH) wavelength is given by |f41| , |42"|| : 

pNL _ pS _|_ -pBQ _|_ pBD ^ 
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where P s is the surface nonlinear polarization, P B Q is the bulk quadruple contribution, and 
P BD is the bulk dipole DCF induced polarization. The last contribution is governed by the 
fourth-rank cubic susceptibility tensor x an d can be written phenomenologically as 



,BD 



X 



{3) ' BD (2u; a, u, 0) : E(cj)E(w)E , (2) 



where E(uj) and Eo are the amplitudes of the fundamental radiation and DCF strength inside 
the semiconductor, respectively For crystals such as Si and Ge of point group symmetry 
Oh, x^' BD has 21 nonzero tensor components f20|l , of which only three are nonequivalent. 



The bulk quadrupole contribution in the plane-wave approximation is given by 

p bq = x m,B Q{2uJ . Wj ^ . E ( w )i k ( w )E( w ) > (3) 

where x^ 2 '' B ® is a fourth-rank tensor which represents the quadrupole contribution to the 
quadratic nonlinear susceptibility from spatial dispersion and k(u;) is the wavevector of the 
fundamental radiation in the semiconductor. x^' E ® has the same symmetry properties as 

X (3),BD_ 

For the surface contribution to the multipole expansion is hardly expected to be 
valid, and we suppose that pE5fl: 



^ S (2 W ; W)W ):E( W )E( W ), (4) 

where x S is a third-rank tensor representing an effective quadratic susceptibility of the 
surface layer, which includes a local part from breaking of inversion symmetry at the surface, 
and a nonlocal part from the discontinuity of the normal electric field component at the sur- 
face. The structure of x S depends on the particular crystalline face under consideration. 

The SH electromagnetic field E(2u>) is found by solving the inhomogeneous wave equation 
for propagation of the SH wave with P NL as a source term f4"5| , |4"l}| . The solution can be 
written formally in terms of the tensorial Green-function G(r, r', 2u), which is defined to be 
the solution of the wave equation with a point source at r'. Since translational symmetry in 
the interface plane is assumed, the SH field is given by 
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E(z, k {{ , 2u) = J G(z, z', fc N , 2to)P NL (z', 2u)dz\ (5) 

where fen is the in-plane component of the SH wavevector. Hereafter we use the xyz coordi- 
nate frame with the xy plane coinciding with the interface and the positive z -axis directed 
toward the semiconductor bulk. Expressions for the components of G, are calculated in 
Refs. [43,44]. The DCF induced part of the SH field is given by: 

E BD (z, 2co) = F 2 ^Xe f %P x 
+00 

x J E (z>) exp (z (k 2 ^ z + 2k Ufg ) z') dz' (6) 


where the scalar factor XePf ls a linear combination of components of x^' BD which depends 
on the experimental geometry, 1^ is the intensity of the fundamental radiation, k W)Z and 
are the normal wavevector components of the fundamental and SHG radiation, respectively, 
in the semiconductor, the unit vector p defines the polarization of the EFISH field, and 
and F 2w are the transmission factors which include Fresnel coefficients and a correction for 
multiple reflection in the silicon oxide at both u and 2uo. Eq.@ properly takes into account 
retardation, the penetration depth of the fundamental wave, the escape length of the SH 
wave and multiple reflection interference effects in oxide layer. 



B. DC-electric-field spatial distribution 

To perform the integration in Eq.(||) one must know the spatial distribution E (z) across 
the SCR. In this section we consider the screening of this DCF within the framework of Fermi 
carrier statistics P5]-^|. The spatial distribution of the electrostatic potential tp(z) in the 



planar semiconductor-dielectric system can be found as a solution of the one-dimensional 
Poisson equation: 

where e is the static dielectric constant of the semiconductor (dielectric) and n = n(z) is the 
space charge density. The boundary conditions for Eq.(|7]) are: 



^(+oo) = /i, 

(p(-D) = fi + (p , (8) 

where /i is the chemical potential of the semiconductor and D is the thickness of the oxide 
film. The first equation in (JD is a statement of charge neutrality in the bulk of the semi- 
conductor. The second equation takes into account the application of external potential ipo 
to the metal electrode with respect to the semiconductor. We divide the charge density into 
field independent and field dependent terms: 

n = n fi + n fd , (9) 

where n/j includes the density of the ionized donors No and acceptors Na, and fixed charge 
n ox trapped in the oxide layer near the semiconductor-dielectric interface: 

n fi = N D + N A + 5(z + 0)n ox , z > 0. (10) 

Hereafter z = +0 and z = — denote positions near the interface just inside the semicon- 
ductor and just inside the dielectric, respectively. 

The spatial distribution nfd{z) is, in principle, a nonlinear functional of the potential 
(p at all points inside the semiconductor. However, first we find expressions for rifd(z) and 
Eq(z) within the model of local screening of the DCF in a Fermi electron-hole gas, in which 
rifd(z) depends on the potential <p at point z, i.e.. nfd((p) = rifd{(p(z)). The field-dependent 
part of charge density consists of the concentration of holes rih, electrons n e , and interface 
traps nu, which depend on the interface potential: 

n fd (z) = n h (<p(z))+n e (<p(z)) + 5(z-0)n it ( V (z = +0)), z > 0. (11) 

Since we assume that the SHG response comes from the semiconductor or semiconductor- 
dielectric interface, we treat charges in the oxide layer as an effective fixed trapped charge 
n ox . Since at z > the variable z does not enter into Eq. (|ICi| , |TTD explicitly and the charge 
density n/d depends on the coordinate via <p(z), the Poisson equation (|7|) has the first 
integral: 
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E 2 M = y / nfoW. (12) 

Using the charge neutrality condition in the bulk of the semiconductor for completely 
ionized donors and acceptors yields 

N D = eiV c $ , (13) 



kT 

Na = -eN v <$> (j^A) ■ (14) 



Eqs.(10,ll) have the form: 



where 



n fl = eN c $ i^r) ~ eN v <S> + 6{z - 0)n ox , (16) 



2 7 

$(r) = -= / v/x (1 + exp (x - r)) 1 dx, (17) 

\/7T 7 



is the Fermi-Dirac integral; Ny and A*^ are the density of states in valence and conduction 
bands, respectively, which depend on density-of-state and effective mass of electrons or holes; 
Ey and Ec are the energies of the upper level of the valence band and the lower level of the 
conduction band, respectively; k is the Boltzmann constant, and T is the temperature. 

Interface traps are charged mid-gap states at the semiconductor-dielectric interface re- 
sulting from interruption of the semiconductor lattice structure or interface imperfections. 
As the interface electrostatic potential changes, the trap levels move up or down while the 
Fermi level remains fixed. Interface trap density is defined in terms of the energy distri- 
bution L A ' D (E) of trap levels across the semiconductor band gap: 

mfo) = ej (L D (E - e(f)F D (^ - E + eip) - L A [E - e<f)F A {fi -E + e<p))dE, (18) 

where superscripts A, D denote acceptor or donor traps, and F A and F D are Fermi distri- 
bution functions: 



F a > d (t)= (l + g A ' D exp (r^y 1 , (19) 

where coefficients g A = 1/4 and g D = 2 reflect the ground-state degeneracy of the ac- 
ceptor and donor levels. The specific form of L A ' D (E) depends on the preparation of the 
semiconductor-dielectric system. In the calculations we model this distribution as a set of 
Lorentz functions. 

Figure 1 shows the distributions (fi(z) and Eq(z) across the SCR of p-doped silicon 
modelled within the above framework. In the depletion regime, where the Schottky approx- 
imation is valid, (p(z) is close to a parabolic function. For larger applied bias corresponding 
to inversion, the SCR divides into a thin subsurface region of rapidly changing potential, 
and a long tail of gradually decreasing potential. The transition depth z between these 
two regions lies at several nanometers. In the accumulation regime, (p(z) drops completely 
within Zq. 



C. The role of surface quantization effects in the subsurface region 

The large gradient of the subsurface <f{z) for accumulation and inversion regimes requires 
that quantum effects be considered in the screening of the DCF. We take them into account 



via self-consistent calculations H48| , using the Hartree-Fock (HF) approach to describe the 
exchange electron interaction. In the following we consider the screening of a "positive" (in 
the above notation) external potential in the subsurface region by electrons, with negligible 
contribution from holes. The opposite case of "negative" potential is treated similarly. 
The HF equation for the single-electron wave function ipi(r) is given by: 

|^(r) + e 2 E < ^(r')|(e|r - r'|)-%(r') >^(r) - 
-e 2 J2 < ^'(r')|(e|r - r'Diifr') >^(r) + U (r)^(r) = E^r), (20) 

jj^i,\\spins 

where Uq(z) = —E z + / ri^(r')(e |r — i*'|) -1 d 3 r' and the sum in the exchange (third) term 
is over states with parallel spins; brackets denote averaging over the stationary state. 
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Because of translation symmetry in the x, y-plane ipi(v) = (p%{z) e ipiT ^. We consider 
the case in which only one energy state for the subsurface electronic -motion is responsible 
for most of the screening. This is confirmed by the numerical results. We also assume 
<Pi(z) = <f(z) to be independent from pi. Then Eq.(pO|) may be written in the form of a 
Schrodinger equation with self-consistent potential U z : 

where 

J e r — r 



n 



r) =< n(r) >= e £ ^(r)V£(r)/(E K ), (23) 



A 



n 



r) is the density operator, f(Ek) — 1/(1 + & ^' kT ) is the Fermi occupation factor, 



e 2 



C(r-r') = l--—— £ / (E,)/(E,)e^^)(-'). (24) 

C(r — r') can be interpreted physically as a correlation function for the in-plane motion of 
electrons. Boundary conditions for the wave function are given by 

£(0) = 0, <p(za) = 0. (25) 

From the equations above one can show that the potential U(z) obeys the following 
equation for the 2D-system under consideration: 



-E + ^ J dz'(n(z')F(z - z') + n fl (z')sgn(z - z')), (26) 



dU „ 2vr 

dz 



where 



= J W{z)C(p\z\) (2?) 
' (1 + p 2 f 2 



The electrostatic potential ip(z) obeys the equation 

^= ~ E ° + yI dz '« z ') + n fi (z'))sgn(z - z'), (28) 
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which can be derived from the Poisson equation. 

Eqs. ( p6] , p8|) differ one from another by the factor F(z) in Eq.(^). It can be shown 
that since C(r) — ► 1 at r — > oo, F(z) —>■ 1 at z — > oo. Therefore, the self- consistent 
potential U(z) is closely related to <f(z). Moreover, remote charge layers contribute equally 
to U(z) and <p(z). Nevertheless, the electrostatic and self-consistent potentials are distin- 
guished by the role of quantum-correlation effects in the electron plasma. The electrostatic 
potential describes interactions of a charged probe particle with other charges only via the 
electromagnetic field. The self-consistent potential U(z) also includes the electron's ten- 
dency to "wedge" itself into other electrons and repulse them via the exchange interaction, 
and therefore differs fundamentally from <p(z). On the other hand, the EFISH bias depen- 
dence is expressed in terms of the classical potential tp(z), because the major contribution 
to the semiconductor optical response comes from bound electrons, whereas screening in the 
semiconductor is caused by free carriers, and there are no correlation effects between these 
two different types of particles. 

Summarizing this section, we have obtained a set of Eqs. ( |2T|J2B|J2^ ) for the self- consistent 
potential U(z), electron wave function ipi( r ) an d charge density n(z). These equations 
describe the screening in the immediate subsurface region z < z . This approach takes into 
account correlation effects in the electronic liquid via the factor F(z) which is related to the 
in-plane correlation function C(r\\ — rj,). 

D. Numerical experiment 

In this section the model bias dependence of the EFISH intensity is found by numerical 
integration of the Poisson equation and the wave equation. First, the dependence of the 
EFISH intensity on the doping of the semiconductor is taken into account. Then, the 
influence of parameters of the semiconductor-insulator interface on the amplitude of the 
EFISH wave is considered. 

To find the DCF induced SH field amplitude E BD for applied bias U, E (z) has been 
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calculated by numerically solving the first integral of the Poisson Equation (0) with the 
charge densities given by Eqs. (|T5|,|ISD . The boundary condition at the surface of the metal 
electrode of MOS structure is given by ipo = U. U is related to the interface field E int = 
E (z = +0) and interface potential ipi nt = ip(z = +0), by 

u = e sc e^ l E int (Lp int )D + cp int . (29) 

The parameters of silicon, which is used as a model semiconductor, have been taken from 
Ref. [49]. According to Eq.(§), the EFISH field E BD is a product of the integral 

+oo 

I(U) =h + il 2 = J E (z) exp (i (k 2 ^ z + 2k Utg ) z) dz } (30) 
o 

and the complex factor F^F^Xeffi which is a bias-independent constant for a given fun- 
damental wavelength. This allows us to neglect the complex value of the latter term and 
simulate the bias dependence of h(U) and hiU) by the bias dependence of ReE BD and 
lmE BD in units of F 2uJ F^Xeff ■ This notation is used in the numerical experiment shown in 
Figs. 2,4 and 5. Figure 2 shows KeE BD and lmE BD as functions of the bias applied to the 
MOS structure, calculated by evaluating the integral in Eq.(^|) for different dopant concen- 
trations of a n-type silicon wafer covered by silicon oxide film 19 nm thick. The fundamental 
radiation wavelength is presumed to be 730 nm. 

Two important trends in these curves are noteworthy. First, lmE BD depends strongly 
on the bias only in the region of negative biases between V and a saturation bias we 
denote as Uq. Outside of this interval the amplitude of the EFISH field saturates. This 
strongly contradicts the previous phenomenological assumption that the amplitude of the 
EFISH field depends linearly (and the EFISH intensity quadrat ically) on the applied bias. 
The saturation of the imaginary part of the EFISH field amplitude for U <Uq and U > is 
attributed to the inversion and accumulation regimes of the external bias screening in the 
SCR (see inset in Fig. 1) as the DCF is mostly localized inside a thin subsurface layer of 
nm-scale thickness. Since the imaginary part of the Green's function is equal to zero exactly 
at the interface, lmE BD becomes practically insensitive to the DCF inside the inversion and 
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accumulation layers. Thus U = Uq and U = define end-points of a bias region which 
corresponds to the depletion regime; the interface potential <po for external bias Uq is equal 
to 2{ei — /i), where is the midgap energy [fjOfl . 

Second, decrease of dopant concentration leads to the decrease of the absolute value of Uq 
and E BD . Figure 3 shows the dependence of the absolute value of Uq on the dopant concen- 
tration of the n-type silicon wafer for various oxide thicknesses. For dopant concentrations 
larger than 10 16 cm~ 3 , the absolute value of Uq scales approximately as the square root of 
No, while for smaller doping levels \Uq\ scales as lniYo, as is clearly shown in the inset of 
Fig. 3. Applied bias Uq according to Eq . (|29|) consists of voltage drops e sc e^ 1 EQ(Lp )D across 
the oxide film, and <y9 across the silicon SCR. Within the Schottky approximation for the 
SCR f4"7| , the interface potential ip and interface field E , corresponding to applied bias Uq 



are given by 

E = 2^& , <pq = 2(6, - kTHNoNc 1 )), £ = 2ireN D e£. (31) 
Therefore, for high doping levels the applied bias mostly drops across the oxide layer and 



Uq(N d ) oc Eq(N d ) oc y Nd ^(NdN^, 1 ). For low doping the interface potential dominates 
and Uq(N d ) oc (po(N D ) oc h^iV^iV,^ 1 ). For thiner oxide layers, less of the applied voltage 
is dropped across the oxide and the transition from a logarithmic to a square root doping 
dependence of Uq occurs at a higher doping level. 

Figure 4 shows the EFISH amplitudes for applied bias Uq as functions of donor con- 
centration, No- Over a wide range of concentrations KeE BD and lmE BD depend on the 
square root of iV^. The latter can be explained by integrating Eq.(^) with a linear DCF 
E(z) = Eq — 2£z across the SCR, as in the Schottky model. This integration yields the 
following expressions for the EFISH field: 



ReE BD oc E A 2 + 2f l , * (32) 



,A| - A 2 
A? + Af 



ImE BD <xE A l -^^^, (33) 
14 



where Ai = Re(2/c w + k^) and A 2 = lm.{2k UJ + fc 2a) ). Since the interface field E = 2£W 
depends linearly on the width of the SCR, the restrictions WA 1 >> 1 and WA 2 >> 1 lead to 



the following expression for the complex SH field: E BD oc -Eo(A 2 + iAi) oc \J No ^{NuNq 1 ). 

Thus ReE BD and lmE BD scale approximately as the square root of N D . Furthermore, 
lmE BD /ReE BD = A X A 2 1 , i.e. the ratio of ReE BD to lmE BD is the ratio of the characteristic 
length scale of absorption, Aj 1 , to that of retardation, A^f 1 , for the SH waves. As the energy 
of the 365 nm SH photon used in our calculations is close to the E 1 critical point of silicon, 
Ai = (1/5.7) nm" 1 and A 2 = (1/21.5) nm -1 are sufficiently large to satisfy the conditions 
WAi » 1 and WA 2 >> 1 for dopant concentrations up to 10 18 cm~ 3 (Fig. 4). As shown 
in the inset of Fig. 4, the ratio lmE BD /ReE BD is close to the value of A : A 2 1 = 3.89. 



A sheet of charged interface states changes the relationship between a potential drop 
across silicon and an applied bias due to the boundary condition for normal components of 
the electric displacement vector D. To demonstrate the role of interface traps in the EFISH 
phenomenon we consider the distribution of trap levels across the silicon band gap as a set 
of Lorentz's functions. The charge density of interface traps n it as a function of the interface 
electrostatic potential is given by 



where M = A, D, j numerates Lorentz functions of the energy distribution of the trap levels. 
Nmj, $M,j and Eqmj denote the effective number of traps per unit area, the width and central 
position of j-th Lorentz peak, respectively. These Lorentz functions simulate the continuous 
energy distribution of traps. By setting Smj — * one can account for discrete levels. 

Figure 5 shows the bias dependence of ReE BD and lmE BD for the MOS structure used 
in our experiment which is comprised of n-type Si with a donor concentration of 10 18 cm" 3 




E. The role of interface states in the EFISH phenomenon 




(34) 



3 
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and a 19-nm thick thermal SiC>2 layer. Interface traps are presumed to be acceptors with 
N A = 10 13 traps ■ cm~ 2 eV~ l and 8a = 0.5 kT. Different central positions are considered: 
£oa = A* (thin lines), Eoa — A* = — 10 kT (thick lines) and Eqa — A* = ~20 kT (dashed lines). 
The distribution of such traps across the silicon band gap is sketched at the inset in Figure 5. 
The dotted lines are presented for comparison to the same field components in the absence 
of traps. For negative biases in the inversion regime, bands are bent in such a way that all 
the trap levels are above the Fermi energy, acceptor traps are empty and bias dependence 
of the SH field components is unaffected by the presence of these uncharged traps. As the 
magnitude of the negative bias is decreased, the bands are bent less and trap levels begin 
to fall below the Fermi energy, first the traps with low energy levels, then those with higher 
energy. Consequently, the bias dependence of ReE BD and lmE BD for low energy acceptor 
traps starts to deviate from the dependence for Na = 0, demonstrating the saturation-like 
feature. This is attributed to the pinning of the Fermi level. As the level of the neutral 
traps crosses the Fermi energy, the charge density of interface traps changes and application 
of a smaller bias leads to a decrease in the voltage drop across the oxide film while the 
interface potential and the DCF spatial distribution remain fixed until the trap level is 
completely filled. The bias dependence of ReE BD and lmE BD for N A ^ passes through 
the SHG intensity zero-point for a flatband voltage Ufb, which depends linearly on the 
interface charge n it . In the case of donor interface traps the same effects are obtained, but 
the bias dependence of ReE BD and lmE BD for N D = and N D ^ differ in the inversion 
regime. 



III. RESULTS AND DISCUSSION 
A. Experimental 

For the EFISH experiments the output of an unamplified Ti-Sapphire laser ranging from 
710 to 800 nm was used. The Ti-sapphire laser generates 120-fs pulses with average power 
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of 200 to 300 mW, which is well below the damage threshold of the semiconductor. The 
p-polarized beam was focused onto the sample at a 45° angle of incidence. Reflected p- 
polarized SHG signal was selected by the use of appropriate filters and directed into a 
photon-counting system. High intensity, high repetition rate, short pulses provided a good 
signal-to-noise ratio in our experiments while avoiding significant sample heating. A small 
split off portion of the fundamental beam was focussed through a z-cut quartz crystal that 
provided a reference SHG signal. 

The MOS structures were fabricated from two types of Si(001) wafers: (I) - a highly- 
doped n-type (10 18 cm -3 , Sb doped ) wafer covered by a 19 nm thick Si0 2 film, and (II) 
- a low-doped p-type (1.5 ■ 10 15 cm~ 3 , B doped) wafer with a 8.7 nm thick Si0 2 film. A 3 
nm semitransparent chromium cap layer, and an ohmic aluminum backside electrode were 
evaporated onto the samples. Single-wavelength ellipsometry was used to measure the Si0 2 
thicknesses. As an independent calibration of the flatband voltage, spatially resolved surface 
photovoltage measurements were performed on the same samples. The external bias voltage 
was applied between the chromium and aluminum electrodes. The SHG response from the 
chromium layer was verified to be negligible in comparison with the SHG signal from the 
buried Si(001)- Si0 2 interface. 

The bias dependence of the rotational azimuthal anisotropy of the EFISH intensity was 
measured over a wide range of the bias voltages at various fundamental wavelengths from 710 
nm to 800 nm. Figure 6 shows the azimuthal dependence of the EFISH intensity measured 
for an n-Si(OOl) MOS structure. The pronounced four-fold symmetric anisotropy of the 
EFISH intensity superimposed on a significant isotropic (i.e. independent from the azimuthal 
angle) background was observed at most biases. Variation of the applied voltage changes 
the amplitudes of both the four-fold symmetric and isotropic contributions, both of which 
increase with increasing the absolute value of the bias. At the center of the applied bias region 
near -2.75 V (upper panel) azimuthal dependence possesses a significant eight-fold symmetric 
component, which appears to be comparable with isotropic and four- fold components for the 
same bias. As the applied voltage passes through this bias the phase of the anisotropy shifts 
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by 7r/4. Similar features of the field- induced rotational anisotropy were observed throughout 
the studied spectral range. Figure 7 shows the azimuthal dependence for a p-MOS structure 
which demonstrates similar behavior, except that the eight-fold symmetric component is 
observed at -1.2 V and the isotropic component is appeared to be quite larger than the 
four- fold one. 

B. EFISH at Si(001)-SiO 2 interface: Role of the spatial DCF distribution 

The azimuthal angular dependence of the SHG intensity from the Si(OOl)- SiC>2 interface 
in the presence of the DCF can be described phenomenologically as optical interference of 
DC-field dependent, isotropic and DC- field independent, four-fold symmetric components of 
the SH field: 

hUi>, V) = \<V) + 6cos(4(^ - Vo))| 2 = 
= c (V) + c 4 (V) cos(4(^ - Vo)) + c 8 cos(8(^ - ^o)), (35) 

where ipo is the azimuthal angle of a maximum of rotational anisotropy, a and b are the 
amplitudes of isotropic and anisotropic components of the SH field. The surface, P s , and 
the bulk DCF induced, P BD , components of the nonlinear polarization, P NL , contribute 
to the isotropic component a while the four-fold symmetric component originates from the 
bulk quadruple polarization, P B Q , For the sake of simplicity we put the amplitude b of the 
four-fold symmetric anisotropic component as a real quantity and define the phase of the 
isotropic component a = a' + ia" with respect to b. As a result the dependence of the EFISH 
intensity on the azimuthal angle, ip, is given by a Fourier expansion (^) with 0-th, 4-th and 
8-th Fourier components: 

Co = a ' 2 + a" 2 + h 2 , c 4 = 2a% c 8 = \b 2 . (36) 

Figure 8 shows the bias dependence of the isotropic Fourier component of the EFISH in- 
tensity (left panel) and of the normalized four-fold Fourier component C4 (2y/2c^j (right 
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panel), which is exactly the a' component of the EFISH field. The eight-fold symmetric com- 
ponent, c 8 , appears to be field- independent throughout the range of applied biases. The 
error bars are the averaged amplitudes of Fourier components c\ and C3 (2y/2csj ■ The 
component Cq(U) is quadratic with a minimum at -3.1 V. The a' component passes through 
zero-point also at about -3.1 V and depends on bias nearly linearly with pronounced de- 
viations from linearity at the edges of the bias range. These bias dependences have been 
fitted within the model described above with the amplitude of the field-independent part of 
a and the flatband voltage Ufb as adjustable parameters. Figure 8 shows the model results 
with Ufb = 0.7 V by solid curves which agree well with the experimental data. The obtained 
value of the flatband voltage significantly differs from either minimum of cq(U) or the bias 
for which a' = 0. This difference is attributed to the optical interference of the DCF depen- 
dent (bulk) and DCF independent (surface) contributions to a. For this highly doped MOS 
structure the entire 8 V range of applied biases of corresponds to the depletion regime. 

Figures 9 and 10 show the bias dependence of c and a' for the p-MOS structure. The 
quadratic behavior of cq(U) with a minimum at -1.25 V and approximately linear dependence 
of a'(U) with deviations at the limits are similar to the trends of the n-MOS structure. 
The model of the experimental data with dopant concentration of 1.5 • 10 15 cm -3 , shows 
a clear step-like feature near the center bias which corresponds to the depletion regime of 
the SCR in p-type silicon. However, such peculiarity has not been observed experimentally. 
This discrepancy between the model and experiment occurs for the surface-quantization 
calculation as well, because at these small biases the surface-quantization effects are not of 
importance. 

One possible explanation for the experimentally measured bias dependences is the influ- 
ence of photoinduced effects on the EFISH intensity. The absorption of femtosecond laser 
pulses leads to the excitation of electron - hole pairs in the SCR. The DCF in the SCR sep- 
arates these photo-induced carriers and the density of the charge injected into SCR for the 
pulse duration r ~ 120 fs is on the order of 10 17 cm~ 3 JJ. The presence of these extra charges 
should lead to a decrease of the SCR width. Systematic theoretical description of this effect 
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requires a model that rigorously accounts for the kinetics of electron-hole recombination 
in the subsurface layer. In our model the parameter responsible for the SCR width is the 
dopant concentration. Therefore, photoinduced effects can be effectively taken into account 
by variation of Np. Thick curves in Figures 9 and 10 show the fit with Np = 10 17 cm -3 , 
which is two orders of magnitude larger than the actual dopant concentration. For such 
a doping level the range of biases which correspond to the depletion regime is sufficiently 
larger (about 4 V) and the transition from depletion to inversion and accumulation occurs 
more gradually than for a lower doping level. Much better agreement of the model with 
experimental data is achieved. 

For biases larger then 4 V, which correspond to the strong inversion regime, clear devi- 
ations of the model from experimental data are obtained. This is attributed to the strong 
localization of DCF inside a very thin subsurface layer where the bulk description of the 
DCF screening is hardly expected to be valid and one should take into account surface 
quantization effects. The dashed curve shows the approximation of the data by the model 
with quantum corrections, which demonstrates a better agreement with experimental data 
points in this bias region. 



C. The EFISH spectroscopy: Bulk origin of DC-field-induced contribution 

Tuning the fundamental wavelength in the vicinity of the direct two-photon E\ transition 
allows measurement of the spectrum of the EFISH intensity and deconvolution of the bulk 
and red-shifted surface contributions to the SHG signal |? |L1 , 51| . Figure 11 shows the bias 



dependence c±(U) for various wavelengths of the fundamental radiation, \ u . Tuning of \ u 
from 800 nm to the two-photon resonance near 3.4 eV (A w = 730 nm) produces stronger 
bias dependence of both c±(U) and cq(U). Further decrease of results in a reduced bias 
dependence. The bulk quadrupole component of the SH field, b = E^j s , contributes to 
both the isotropic, cq, and the four-fold symmetric, C4, Fourier components. To extract the 



spectral dependence of the EFISH field, E BD , one must find the spectrum of 



F BQ 

anis 



. The 



20 



latter has been obtained from the spectrum of the eighth Fourier component cs averaged 
over the entire bias region. Integration of the product of the Green's function and the bulk 



quadrupole polarization, according to Eq. 
form: 



gives the spectral behavior of 



E, 



BQ 



X 



m,BQ 



(n) 



in the 



(37) 



Hereafter, x^' B ®(ty is the magnitude of a combination of ^( 2 )> B( 2 tensor components re- 
sponsible for the four-fold symmetric part of P B< 2. J w is the fundamental intensity. Figure 12 



X 



(2),BQ 



(O) 



shows the spectrum of the magnitude of the effective quadruple susceptibility 
The filled symbols in Figure 12 show the spectral dependence of the effective cubic sus- 
ceptibility ^( 3 )' B( 3 extracted from the set of the bias dependences cq(V, A). Both spectral 



dependences of 



X 



(2),BQ 



(O) 



and 



X^' BD {&) peak at approximately 3.4 eV and have been 



fitted by a single Lorentz function with a real spectral background: 

P 



x M (n) = a + 



(38) 



X 



(2),BQ 



(n) 



and 



Q — ujm + iS 

with M = BQ, BD. The solid curves in Figure 12 show the spectral fits of 
t^( 3 )> B£) (0) by Eq. (|38D with the parameters presented at Table I. The values of resonance 
positions obtained are shown to be close to 3.38 eV. This is consistent with the energy of the 
bulk Ei critical point as known from linear spectroscopy and fully indicates a bulk origin of 
the EFISH response. 



D. Low-frequency electromodulation SHG spectroscopy of Si(001)-SiO2 interface 



Modulation techniques are widely used in optical spectroscopy |52| because of their sensi- 
tivity. The right side of Fig. 13 shows the schematic of the low-frequency electromodulation 
of the SHG signal from Si-Si02 interface in a MOS structure by the application of the super- 
position of DC-bias U and low-frequency squarewave modulation voltage AU(Q) with the 
amplitude AU and frequency Q. Microwave frequency and pulse-voltage modulation of the 
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SHG response in Si-based MOS structures were studied in Refs.[37] and [33], respectively. 
Low-frequency electromodulation SHG from GaN surface in electrochemical cell was studied 
in Ref. [31]. 

The efficiency of the modulated SHG signal a{U, Vt) at certain DC-bias U can be de- 
fined by a relative increment of the EFISH intensity while applying the modulation voltage 
AU (VI): 

n(TT , 2(I 2 ^(U + AU)-I 2uJ (U-AU)) dI 2uJ (U) AUjO) 

° {U ' ^ = (I 2 .(U + AU) + I 2 ^U-AU)) w ~dU (39) 

and appears to be a differential characteristic of the EFISH phenomenon which is comple- 
mentary to the static EFISH dependence I 2u> (U). 

Figure 13 shows the experimental static (DCF-induced) EFISH bias dependence 
hu (U + AU) measured at p-Si(OOl) MOS structure for AU = 0.6 V which is typically 
featureless for the DCF induced SHG. Figure 14 shows the experimental bias dependence 
a(U) (open symbols) of the efficiency of the modulated EFISH for for AU = 0.6 V, and 
Q = 100 Hz at an azimuthal angle ip = that minimizes the anisotropic EFISH intensity. 
This dependence shows a flat feature in the vicinity of U = V. This feature is not seen on 
the bias dependence of the numerical derivative of I 2u> (U) function shown in Fig. 14 (solid 
symbols). The bias dependence a(U) as calculated in the framework of the phenomeno- 
logical model above (Sections A and B in Part II) is presented in Fig. 14 (solid curve). 
This curve does not show a flat feature in the vicinity of U = 0. This flat-like feature in 
the experimental differential bias dependence is likely coming from the large modulation 
amplitude, which could not be reduced for technical reasons. 

The inset in Fig. 12 shows the spectral dependence of the modulation efficiency in the 
tuning region of the Ti: Sapphire laser. A peak in the spectral dependence a(X CJ ) is observed 
at the two-photon energy 2hu = 3.41 eV with half-width hAu = 0.023 eV. The spectral 
position of this peak is close to the bulk E 1 resonance. This confirms once again the bulk 
origin of the DCF-induced term of the nonlinear polarization in Eqs. (1) and (2). The 
spectral half-width of the resonance in the differential response a^X^) is smaller then the 
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half-width of the resonances of the electrostatic (DCF-induced) EFISH terms in Fig. 12. 
This shows the increased sensitivity of the EFISH modulation spectra to density of states 
in the semiconductor valence and conduction bands. 

IV. CONCLUSIONS 

In summary, the DC-electric-field-induced SHG and the low-frequency electromodulation 
SHG spectroscopy of Si(001)-SiO 2 interfaces in p- and n-type Si(001)-SiO 2 -Cr MOS struc- 
tures have been studied. The dependences of the DC-electric-field-induced SHG intensity 
on the applied bias are shown to be sensitive to the doping concentration of silicon, oxide 
thickness and fundamental and SHG wavelengths. From spectroscopy of the anisotropic 
EFISH dependences the field-induced contribution has been extracted and the spectrum of 
the cubic susceptibility appears to be peaked at the energy of the bulk E\ critical point. 
The presence of significant EFISH contribution at an unbiased Si(001)-SiO2 interface due to 
the initial band bending has been observed. This initial band bending contribution should 
be taken into account in the further interpretation of the spectroscopic SHG measurements 
at Si(001)-SiO 2 interfaces |,||. 

A general phenomenological model of the EFISH phenomenon is developed. This includes 
a comprehensive analyses of the generation and the propagation of the EFISH wave in the 
silicon space charge region taking into consideration the retardation and absorption effects, 
optical interference of the DC-field-dependent and DC-field-independent contributions to the 
SH waves and interference of multiple reflections in the oxide layer. The spatial distribution 
of the DC-field- induced bulk dipole nonlinear polarization is calculated using the rigorous 
DCF distribution across the SCR taking into account surface quantization effects. The 
influence of the silicon doping level, oxide thickness, interface states and oxide charge traps 
on the screening of the external DCF in the SCR is studied. We have demonstrated the 
sensitivity of the EFISH probe to the charge characteristics of the Si(001)-SiO2 interface 
which makes this technique promising as a noninvasive sensor of the MOS devices for the 
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mapping of interface charge distribution. 
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Fig. 1: The spatial electrostatic potential (left panel) and DC-electric field (right panel) 
distribution across the SCR of p-doped silicon (doping concentration of 1.5 • 10 15 cm~ 3 ) for 
different values of interface potential: +0.95 V (inversion), +0.6 V (depletion) and -0.33 V 
(accumulation). The upper panel is the sketch of the potential and field distribution across 
the MOS structure. 

Fig. 2: The bias dependences of the real (left panel) and imaginary (right panel) parts 
of the EFISH field E BD for different doping levels of silicon wafer. Parameters of MOS 
structure used are in the text. The flatband voltage is supposed to be zero. Insets: the bias 
dependences in the vicinity of zero-point bias. 

Fig. 3: The absolute value of the depletion bias U vs. the doping concentration N D of the 
n-Si wafer for Si- SiC>2- metal MOS structures with different Si02 thicknesses: 1 nm (filled 
squares), 8.7 nm (open squares), 19 nm (filled circles), and 50 nm (open circles). Solid 
curves are guides to the eye. Inset: dependence U (N D ) for MOS structure with 1 nm thick 
oxide in the linear scale. 

Fig. 4: The real (open circles) and imaginary (filled circles) part of E BD for the depletion bias 
Uo vs. the doping concentration No of the silicon wafer with 19 nm thick oxide film. Solid 
curves are guides to the eye. Inset: the doping dependence of the ratio of lm(E BD ) /Re(E BD ) . 
Fig. 5: The bias dependences of the real (left panel) and imaginary (right panel) part 
of the EFISH field E BD for different parameters of acceptor interface states. The en- 
ergy spectrum of interface states is simulated by Lorenzian function with density of 
Na = 10 13 traps ■ cm _2 eV _1 , width of 5a = 0.5 kT, and different central positions - Sqa = ^ 
(thin curves), e a — A* — —10 kT (thick curves) and e a — A* = — 20/cT (dashed curves) 
and sketches on the inset. The dotted curves present for comparison the same dependences 
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without traps. Parameters of the MOS structure are in the text. 

Fig. 6: p-in,p-out SHG signal from n-Si(OOl) MOS structure at several biases for A^ = 
725 nm (2hu> = 3.43 eV) vs. sample azimuthal angle. Solid curves are fits to data by the 
0-th, 4-th and 8-th Fourier components. 

Fig. 7: p-in,p-out SHG signal from p-Si(OOl) MOS structure at several biases for for A w = 
730 nm (2hcu = 3.41 eV) as a function of sample azimuthal angle. Solid curves are fits to 
data by the 0-th, 4-th and 8-th Fourier components. 



Fig. 8: Bias dependences of isotropic cq and normalized four- fold a' = c 4 (2y/2^) SHG 
Fourier amplitudes from n-Si(OOl) MOS for for A w = 725 nm (2fUv = 3.43 eV). Solid curves 
are fits to data using the model presented. 

Fig. 9: Fig. 8. The isotropic SHG component from p-Si(OOl) MOS structure for A w = 730 nm 
(2hu = 3.41 eV) as a function of applied bias. Curves are fits to data using the model of 
the DCF screening within "classical" approach for Na = 1.5 • 10 15 cm~ 3 (thin curve) and 
N A = 10 17 cm~ 3 (thick curve) and with surface quantization corrections (dashed curve). 



Fig. 10: The normalized four-fold symmetric SHG component a' = c 4 [2y/2c$ j from p- 
Si(001) MOS structure for for A^ = 730 nm (2hw = 3.41 eV) as a function of applied bias. 
Curves are fits to data using the model of the DCF screening within "classical" approach 
for Na = 1.5 • 10 15 cm -3 (thin curve) and Na = 10 17 cm~ 3 (thick curve) and with surface 
quantization corrections (dashed curve). Inset: voltage dependences near zero-point of bias. 
Fig. 11: The bias dependences of the four-fold symmetric anisotropic SHG component c 4 (£7) 
for several wavelengths of the fundamental radiation and their fit presented by solid lines. 
Fig. 12: The spectral dependence of module of cubic dipole and quadratic quadruple sus- 
ceptibility around the direct two-photon E 1 transition extracted from the spectra of the 
EFISH azimuthal dependences. Solid lines are fits to data by the Lorenz function with 
real background. Inset: the spectral dependence of the efficiency of modulated EFISH for 
AU = 0.6 V, n = 100 Hz in p-Si(OOl) MOS structure. 

Fig. 13: The experimental static bias dependence of the SHG intensity. The right side shows 
the schematic of the low-frequency electromodulation of the SHG signal by the application 





26 



of the superposition of DC-bias U and low-frequency squarewave modulation voltage AU (fl) 
with the amplitude AU and frequency Q. 

Fig. 14: The bias dependence of the efficiency of modulated EFISH signal in p-Si(OOl) MOS 
structure: open symbols are experimental dependence of a(U) for for AU = 0.6 V, f2 = 
100 Hz; solid symbols are numerical derivative of the static (DCF-induced) bias dependence 
of the EFISH intensity; solid line is the model calculation in accordance of Sections A and 
B in Part II. 
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Z(nm) 




-0.5 H 1 1 1 1 1 1 1 1 1 1 1 

-20 -10 -10 10 

Applie d Bias (V) 




TT| 1 1 I I | 1 1 I I | 1 1 I I | 1 1 I I | 1 1 I I | 1 1 I I | 1 1 I I | 

10 12 10 13 10 14 10 15 10 16 10 17 10 18 10 19 

Doping Concentration (c m" 3 ) 




1 I 1 I 1 I 1 ' 1 ' 1 ' 1 

-20 -10 -10 10 

Applied Bias (V) 



Azimuthal Angle ( d e g .) 

90 ISO 270 360 




Azimuthal Angle (deg.) 



Azimuthal Angle ( d e g .) 

90 ISO 270 360 

_l i I i I i I i L 




90 180 270 90 180 270 360 

Azimuthal Angle (deg.) 




Applied Bias (V) 



Q=100 Hz 
A£/=0.6 V 




U(t)=U +AU(Q) 



I L. 



2 2 4 6 
Applied Bias (V) 
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